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ABSTRACT 

We first compare the CMB lensing model of Seljak (1996) with the empirical model 
of Lieu & Mittaz (2005) to determine if the latter approach implies a larger effect on 
the CMB power-spectrum. We find that the empirical model gives significantly higher 
results for the magnification dispersion, a, at small scales (9 < 30') than that of Sel- 
jak (1996), assuming standard cosmological parameters. However, when the empirical 
foreground model of Lieu & Mittaz is modelled via correlation functions and used in 
^ . the Seljak formalism, the agreement is considerably improved at small scales. Thus 

' we conclude that the main difference between these results may be in the different as- 

£T) \ sumed foreground mass distributions. We then show that a foreground mass clustering 

CO ■ with £(r) oc r -3 gives an rms lensing magnification which is approximately constant 

0^ ' with angle, 8. In Seljak's formalism, we show this can lead to a smoothing of the CMB 

power-spectrum which is proportional to e~ a 1 I" 1 and which may be able to move the 
first CMB acoustic peak to smaller I, if the mass clustering amplitude is high enough. 
Evidence for a high amplitude of mass clustering comes from the QSO magnification 
& • results of Myers et al (2003, 2005) who suggest that foreground galaxy groups may also 

Oh' be more massive than expected, implying that Q m ~ 1 and that there is strong galaxy 

anti-bias, b « 0.2. We combine the above results to demonstrate the potential effect 
' of lensing on the CMB power spectrum by showing that an inflationary model with 

^ \ neither CDM nor a cosmological constant and that predicts a primordial first peak 

at I = 330 might then fit at least the first acoustic peak of the WMAP data. This 
model may be regarded as somewhat contrived since the fit also requires high redshift 
reionisation at the upper limit of what is allowed by the WMAP polarisation results. 
But given the finely-tuned nature of the standard ACDM model, the contrivance may 
be small in comparison and certainly the effect of lensing and other foregrounds may 
still have a considerable influence on the cosmological interpretation of the CMB. 
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1 INTRODUCTION 

The standard ACDM model gives an excellent fit to the 
CMB power spectrum from WMAP and other experiments 
(Hinshaw et al. 2003, 2006). The reduced chi-squared val- 
ues are impressive for a model with 7 basic parameters and 
fitting over 40 approximately independent data points. The 
model also appears to give an excellent fit to independent 
datasets such as the 2dFGRS galaxy power spectrum (Cole 
et al. 2005). 

Impressive though these fits are, the standard model 
still appears to have a variety of more fundamental prob- 
lems (see Shanks, 2005 and references therein). For exam- 
ple, the CDM particle still has its own associated problems. 
The limits from the Cryogenic Dark Matter Search (Akerib 



et al, 2004) and other experimental constraints when com- 
bined with the WMAP result for the CDM particle density, 
now leave only a small amount of 'generic' parameter space 
allowed by some of the basic supersymmetric models for the 
neutralino (e.g. Ellis et al 2003, Barenboim & Lykken 2006); 
this is known as the 'neutralino overdensity' problem. There 
is also what is known as the 'gravitino problem'. If the neu- 
tralino is the lightest supersymmetric particle and the grav- 
itino is the next-to-lightest then if the neutralino is stable 
as required for a CDM candidate then the gravitino may 
be also massive and unstable; if its lifetime is longer than 
0.01 sees then the decay of these particles into other par- 
ticles such as baryons could upset the Big Bang prediction 
of the light element abundances (e.g. Cyburt et al 2003). 
Thus the usual agreement of light element abundances with 
the baryon density may not be taken for granted if the neu- 
tralino is the CDM particle. In some sense, the invoking of 
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the CDM particle could appear to have over-specified the 
model in terms of solving the baryon nucleosynthesis prob- 
lem. 

The cosmological constant or dark energy presents fur- 
ther fundamental problems. Indeed, string theory prefers a 
negative A rather than the observed positive A (e.g. Witten, 
2001). But it is the small size of the cosmological constant 
that undoubtedly presents the biggest difficulty for anyone 
assessing how successful the standard model is in explaining 
the current data. In the standard model, just after inflation 
the ratio of the energy density in radiation to the energy 
density in the vacuum was one part in 10 120 . This means 
that the standard model ends up potentially having more 
fine-tuning than implied in the original 'flatness' problem 
before it was solved by inflation (Guth 1981). It is this prob- 
lem that presents the difficulty in assessing the success of the 
standard model where the low reduced chi-squared, from fits 
to the CMB power-spectrum has to be balanced against the 
huge degree of freedom allowed by the invocation of a very 
small cosmological constant. 

Here we shall take the view that the fundamental com- 
plications of the present standard model suggest that it may 
still be worthwhile to explore if other models may explain 
the data. This may imply the introduction of further model 
parameters which may locally increase the model's compli- 
cation, but if they then allow the dropping of the hypothesis 
of the cosmological constant and even the CDM particle, the 
model's complication will be globally decreased. 

CMB foregrounds will form our main escape route from 
the tight constraints set by the CMB power-spectrum. One 
example will be the evidence for the effects of high red- 
shift (z « 10) reionisation from the WMAP polarisation 
results (Kogut et al 2003, Page et al 2006). This already de- 
creases the height of the acoustic peaks by £s 20% relative 
to the multipoles at larger scales even if the reionisation is 
homogeneous. But galaxy clusters in the CMB foreground 
at lower redshifts provide several other ways of contami- 
nating the CMB signal. The hot gas in the richer clusters 
scatter the CMB photons via the SZ effect and there are 
claims from the cross-correlation analysis of WMAP data 
with Abell clusters that the SZ decrements may extend up 
to 0.5 degree scales (Myers et al 2004) . They can also grav- 
itationally affect the CMB photons primarily via the Sachs- 
Wolfe effect and also by gravitationally lensing the CMB 
photons. Although the effects of foreground lensing are gen- 
erally claimed to be small, we shall start by considering the 
claims by Lieu & Mittaz (2005) that the effects of lensing 
on the CMB power spectra may be bigger than previously 
estimated. Lieu & Mittaz analysed lensing by foreground 
non-linear structures by using catalogued properties of clus- 
ters and groups rather than power-spectra, modelling the 
clusters as singular isothermal speheres. 

We shall then relate the above approach to the CMB 
lensing formalism of Seljak (1996) as implemented in CMB- 
FAST (Seljak & Zaldarriaga 1996). We shall show that the 
standard formalism also allows models with a first peak at 
larger wavenumber, I, than the standard model to fit the 
WMAP data, if the mass power spectrum has a strong anti- 
bias with respect to the galaxy power spectrum as suggested 
by the QSO lensing results of Myers et al (2003, 2005). In 
this way we shall show that a large range of models may be 
able to reproduce the observed first peak. In particular, wc 



shall argue that a low Ho, Qbaryon = 1 model may give an 
excellent fit to the CMB TT power spectrum with no need 
to invoke either a CDM particle or dark energy. 

Hence, in Section 2 we review and modify the CMB 
lensing results of Lieu & Mittaz. In Section 3 we introduce 
the CMBFAST lensing formalism and compare with the re- 
sults of Lieu & Mittaz. In Section 4 we use these results 
together with cosmological parameters determined from the 
QSO lensing results of Myers et al (2003, 2005) to predict 
a significantly increased effect of lensing on the CMB. In 
Section 5 we then consider the baryonic cosmological model 
that has an intrinsic CMB first peak at I w 330. We discuss 
how homogeneous reionisation at high redshift and then the 
effect of lensing can significantly improve agreement with 
the observed first peak at I w 220. In Section 6 we discuss 
our conclusions. 



2 LENSING MODEL OF LIEU AND MITTAZ 
(2005) 

We shall first consider the lensing model of Lieu & Mit- 
taz (2005). These authors have suggested that the effects 
of lensing by foreground groups and clusters on the CMB 
can be much larger than previously calculated. This area 
has been controversial in the past with some authors claim- 
ing large effects of lensing on the CMB (Fukushige et al., 
1994, Ellis et al., 1998) while others find much smaller ef- 
fects (eg Bartelmann & Schneider, 2001). The lensing ef- 
fects claimed by Lieu & Mittaz are significantly larger than 
those claimed by Bartelmann & Schneider and somewhat 
larger than claimed by Seljak. This may be because Lieu 
& Mittaz use strong-lensing formulae appropriate for clus- 
ters modelled as isothermal spheres or NFW profiles rather 
than the weak-lensing approximations used by Seljak and 
Bartelmann & Schneider, although all the above authors 
assume the weak-lensing Born approximation in terms of 
considering only small deviations from the original photon 
path. Another reason may be that they use empirically de- 
termined models for the foreground structure which may 
produce larger effects than standard ACDM models. 

Lieu & Mittaz treat the demagnification in the empty 
voids and the magnification in the clusters separately and 
show that in terms of the average magnification (rj) they 
cancel, so that the same 'light conservation' result applies in 
this inhomogeneous case as in the homogeneous case treated 
e.g. by Weinberg (1976). They show that the same result 
applies when strong lensing is taken into account. Here we 
follow Lieu & Mittaz and define r\ = (6' — 6) /0 as the magni- 
fication 'contrast'; the effect of the scattering is to increase 
the average angular size 6 of the source by this fractional 
amount. 

The basic relations are for the expected value of rj in the 
case where the source is far behind the lens from equation 
(18) of Lieu & Mittaz:- 

(ti) = |n groupSJ ff 2 / ' dx [i + z(x)] {x " ~ x)x , (1) 

z Jo Xs 

and equation (31) of Lieu & Mittaz gives for Sr/ 

(Svf = 5f no*? (l - g + g) J] W ^,U)(2) 
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where 

/(^,6 min ) = ln(-^)- A (3) 

V O m in / 7T^ 

and Pij is the probability of finding a group with velocity 
dispersion <t; and radius Rj . Q 

In one case, they use a space density of groups which 
corresponds to that observed in the group sample of Ramella 
et al (1999, 2002) and assuming the other group parameters 
such as the velocity dispersion, a. This produces (rj) « 0.099 
and Sr/ w 0.093, assuming the groups continue with the same 
comoving space density out to Zf = 1 in a fi m =0.27 model. 
Lieu & Mittaz note that this assumption may not be justi- 
fiable for groups because of the effects of dynamical evolu- 
tion. For their Abell cluster model they find (rj) m 0.099 and 
Sr) — 0.064. Lieu & Mittaz argue that their assumption of no 
dynamical evolution of the cluster density may have more 
empirical support in the case of clusters than groups. One 
point to note is that Lieu & Mittaz use the line-of-sight rms 
velocity dispersion and the full 3-D rms velocity dispersion 
should be \/3 bigger. Although (r/) would remain unchanged, 
Sr) would then increase by a factor of 3, implying Sr/ = 0.28 
for the Ramella et al groups and Sir = 0.19 for the Abell 
clusters. 

These values for Sr/ apply in the coherent regime of scat- 
tering where a CMB patch has a mean angular size less 
than a, the size of a group at an average redshift. At larger 
scales, in the incoherent regime of scattering, the effective 
Sr] reduces according to Srj oc 1/8 (see equations 33, 34 of 
Lieu & Mittaz). For the groups and clusters considered by 
Lieu & Mittaz, the division between these regimes occurs at 
m lOarcmin. The form of Sr] (= a(9)/9 in the terminology 
of Seljak, 1996) with 6 is given in Fig. 1 for the above mod- 
ified group and cluster models of Lieu & Mittaz. We shall 
compare these to the standard CMBFAST estimates of the 
magnification function in the next section. 

Lieu & Mittaz note that in cases where Sr) is large com- 
pared to (rj) a skewed distribution for r) may be implied due 
to negative density fluctuations cutting off at p = (Dyer- 
Roeder, 1972, empty beam). Although this may also apply 
to the models above and elsewhere in this paper, here we 
shall simply assume a symmetric magnification distribution 
which is applied using equation A7 of Seljak (1996). 



3 LENSING MODEL OF SELJAK (1996) 

We next discuss the more usual lensing formalism of Sel- 
jak (1996) as implemented in CMBFAST. This is based on 
inferring the fractional magnification = <j(8)/8 (= Sr} in 
the terminology of Lieu & Mittaz) from the mass power 
spectrum using the variation of Limber's formula given in 
equations (7) of Seljak (1996). Using the CMBFAST source 
code, a (6)/ 6 can be generated. The ACDM model with 
Q m — 0.27, £Ia = 0.73 (as assumed by Lieu & Mittaz) with 
the non- linear extension of Peacock & Dodds (1994) gives 
the result shown in Fig. 1. It can be seen that the a(ff)/9 
drops like 1/9 at large scales as expected in the incoherent 
regime while flattening at smaller scales. At 9 > lOOarcmin, 

1 Note that equation (32) of Lieu & Mittaz (2005) contains a 
typographical error; the numerical value should read 2.7 X 10 . 
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Figure 1. Lensing dispersions, 8rj = a/6, for the models of Lieu 
& Mittaz and Seljak (1996). The LM clusters and groups mod- 
els are those of Lieu & Mittaz where the lenses are respectively 
Abell clusters and the groups of Ramella et al in a ACDM cos- 
mology (see Sect. 2). The results of Lieu & Mittaz have been 
multiplied by 3x to base the results on the full 3-D rms velocity 
dispersions rather than line-of-sight dispersions. The LM groups 
(Myers) model assumes each group contains the high mass im- 
plied by the QSO lensing results of Myers et al. and assumes an 
EdS cosmology (see Sect. 4.1). The ACDM and SCDM models are 
the results from CMBFAST including the non-linear extensions 
from Peacock & Dodds (1994) and are comparable to the results 
of Seljak (1996). The models based on the power-law correlation 
functions also use Seljak's formalism; the power-law shown with 
index -3 assumes ro = 7.4h _1 Mpc and the -1.8 power-law assumes 
ro = 6. 0h — x Mpc. The constant a = (=0.005 rad.) model has an 
amplitude similar to that implied by the QSO lensing results of 
Myers et al. 

the ACDM model lies between the group and cluster results 
of Lieu & Mittaz. However, at smaller scales, the models of 
Lieu & Mittaz are 2 — 3x higher than the ACDM model. 
Tests suggest that the neglect of linear theory evolution in 
the models of Lieu & Mittaz does not appear to explain this 
discrepancy. An SCDM model is also shown in Fig. 1. With 
COBE normalisation, this model gives erg = 1.6 and lies 
above the ACDM model at all scales and above the models 
of Lieu & Mittaz at large scales. 

The magnification rms dispersion, <r, is then used in 
Seljak's equation A6 which, under the assumption that the 
isotropic term is dominant, becomes his equation A7: 

/*oo 

C(9) = (2ty)- 1 ldle-" 2mi2/2 CiJ (l9). (4) 
Jo 

where C; is the CMB angular power spectrum and C(ff) is 
the lensed angular correlation function. This equation has 
the appearance of a Hankel transform but the dependence 
of o on 9 means that this is only true under the assumption 
that the 9 dependence is slow. 

When equation A6 is used with the standard model, 
only small effects are seen on the CMB (see Seljak, Fig. 2). 
Indeed, we found that any a(9)/8 that flattens at scales 
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smaller than 1 degree even with amplitudes as high as 
a(9)/9 — 0.45, generally gave small lensing effects on the 
first acoustic peak. Since the lensing dispersion is close to 
that of Lieu & Mittaz, similar results are found under their 
assumptions. 

A case of interest is the one where we start by assum- 
ing that a(6) is a constant. This means that a(9)/9 oc 1/9 
i.e. the usual large-scale behaviour extends to the smallest 
scales. It also means that the combination e _CT ^ l ^ 2 Ci is 
the lensed power spectrum. By adjusting the size of a it is 
then possible to smooth small scale peaks away and even 
to partially smooth away and hence move the first peak. 
For the moment, we simply assume a = 0.005 rad. in the 
case shown in Fig. 1 for 9 > 7'. The question is whether 
this amplitude and form of a(9)/9 can be justified physi- 
cally. Equation 7 of Seljak(1996) implies that a =constant 
requires the power spectrum of the density, P{k), (oc P^, x k 4 
where is the power spectrum of the potential) to have the 
form P(k) ~ k in the range of interest which implies that 
£(r) ~ r~ 3 , since Fourier transforming the density power 
spectrum, P(k) oc k n , to the correlation function implies 
£(r) oc r ~ (3+n) asymptotically for -3 < n < 0. Thus it ap- 
pears that the constant a model will require a correlation 
function for the mass which is steeper than the observed 
galaxy correlation function, £(r) oc r~ 1,8 . As can be seen, 
the amplitude of a is also much higher than for the stan- 
dard model. But we shall see next that such an amplitude 
may be suggested by considerations of QSO lensing and that 
the a ^constant form shown in Fig. 1 then can be used to 
'move' the first acoustic peak. 

We now further draw the comparison with the results of 
Lieu & Mittaz by expressing their model in power-spectrum 
or correlation function terms. These authors assume a ran- 
dom distribution of Abell clusters each with a singular 
isothermal sphere for its mass distribution. According to 
Peebles (1974) a model with randomly placed monolithic 
clusters with density profile p(r) oc r~ e implies a correla- 
tion function with £(r) oc r~( 2e ~ 3 K e = 2 as assumed by 
Lieu & Mittaz gives 7 = 1.0 whereas e = 2.4 gives 7 = 1.8. 
Fourier transforming the correlation function to the power 
spectrum (as above) we find that isothermal spheres of large 
radius leads to a power spectrum of slope k~ 2 which can be 
assumed to apply in the range 0.01 < k < 10 h Mpc -1 . 

When we take the value for the local galaxy correlation 
function scale length of ro = 6h _1 Mpc we find that, assum- 
ing e = 2.4 for the Abell cluster model, a(ff)/9 is approxi- 
mately the same at large scales as the linear theory A model 
(see Fig. 1). Thus, as expected, in the Seljak formalism an 
approximately unbiased mass correlation function with a - 
1.8 power-law gives approximately the same result for a/9 
as the ACDM model with non-linear extension. 

If we now assume a model in the spirit of Lieu & Mit- 
taz where all the mass is in galaxies and all the galaxies 
are in randomly distributed Abell clusters, this will give a 
mass correlation function that extends as a -1.8 power-law to 
> 10h _1 Mpc with a correlation length of ro = 10.6h~ 1 Mpc 
(from equation 23 of Peebles, 1974). This then produces 
a a/9 result which would be 1.7x higher than that for 
ro = 6h _1 Mpc shown in Fig. 1. This model gives a gen- 
erally flatter result with 9 than the LM cluster (or group) 
model, agreeing with the LM cluster model at 1 arcmin, ly- 
ing a few times lower at 10 arcmin and a few times higher 



at 100 arcmin. The cut in the cluster mass profile at radius 
w 2h _1 Mpc assumed by Lieu & Mittaz makes the compari- 
son less valid at larger scales. An w 2.5 x bigger correlation 
function scale length would be required to reach the con- 
stant a model even at large scales ie ro w 27h~ 1 Mpc. Also, 
as we have seen the model implied by the -1.8 power law for 
£ is significantly flatter than the constant a model discussed 
above and the SIS model with e = 2 and 7 = 1 as actually 
assumed by Lieu & Mittaz would give an even flatter slope. 

We conclude that Lieu & Mittaz and Seljak appear 
to obtain reasonably similar results for a/9 when simi- 
lar input models are assumed. Although the Lieu & Mit- 
taz group/cluster models show more power at small scales 
than the ACDM model, when we assume mass correlation 
functions which should be reasonably appropriate for the 
Lieu & Mittaz cluster model in the Seljak formalism, im- 
proved agreement is seen at small scales. Therefore in terms 
of the discussion (e.g. Lewis & Challinor 2006) of how the 
results of Lieu & Mittaz relate to earlier work, there may be 
no serious disagreement, with their assumed higher ampli- 
tude of foreground mass clustering probably being a bigger 
factor in explaining any difference than their partial inclu- 
sion of strong lensing in their lensing formalism. But in the 
wider context of affecting the first peak position, such dif- 
ferences between the models are academic because none of 
the a(9)/9 models which flatten at small scales produces 
the simple smoothing of the CMB first peak of the constant 
a model; models with more small-scale power and a larger 
mass clustering amplitude are required. We next look to see 
if there is any evidence that the size of the magnification 
dispersion may be currently underestimated. 



4 QSO MAGNIFICATION RESULTS 

4.1 Group masses from QSO lensing applied to 
results of Lieu & Mittaz. 

We now consider the possible implications if galaxy group 
masses were as big as those found in the QSO magnifica- 
tion results of Myers et al (2003, 2005). Myers et al (2003) 
found that the QSO-group cross-correlation function fitted 
SIS mass profiles with a — 1125kms _1 . This velocity dis- 
persion is comparable to the values found for rich clusters, 
although the sky density of these groups is w ldeg~ 2 com- 
pared to « 0.1deg~ 2 for the Abell clusters. Myers et al 
(2005) also found that 2QZ QSO lensing by individual galax- 
ies again produced more anti-correlation than expected from 
the standard model. Here they analysed the results using the 
galaxy power spectrum after Gaztanaga (2003) and found 
that in the standard cosmology an anti-bias with b w 0.1 on 
« lh _1 Mpc scales would be required to explain the strength 
of anti-correlation seen. Guimaraes, Myers & Shanks (2005) 
showed that these results were not due to selection effects 
by showing, for example, that the same cluster finding algo- 
rithms in the Hubble volume N-body simulation produced 
much smaller lensing effects. 

Now we discuss below other observations which disagree 
with the results of Myers et al (2003, 2005). Nevertheless, 
here we consider the possibility of whether galaxy groups 
may be dynamically young and the virialised assumption 
therefore unjustified. In this case the group lensing masses 
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may be more correct, implying that the group masses are 
bigger than expected and comparable to the mass of Abell 
clusters. If the full factor of ten increase observed by Myers 
et al for QSO lensing over what is expected for the concor- 
dance model is translated directly into the results of Section 
2 for CMB lensing, then w 3x bigger values of Srj will be 
found, resulting in <5n w 0.6 for both groups and clusters in 
the above model. 

However, Myers et al pointed out that their space den- 
sity of high mass groups implies a Universe with a mass 
density more consistent with that of an Einstein-de Sitter 
model, rather than the ACDM model. The space density of 
groups reported by Myers et al is no = 3±lx 10" 4 h 3 Mpc -3 , 
close to the value of Ramella et al as used by Lieu & Mit- 
taz, (n = 4.4 x 10~ 4 h 3 Mpc~ 3 ). Taking a radius of2h _1 Mpc 
with 6 m in = 7h~ kpc and substituting these in the SIS re- 
sult for <5n, we find 5rj = 0.15. Here we have assumed that 
the group and cluster distribution is now unevolved out to 
Zf = 0.5; our cluster and group densities are essentially de- 
termined within this z range, reducing uncertainties due to 
evolution. We note that Lieu & Mittaz find an average value 
of the group velocity dispersion of 270kms~ 1 whereas the 
weighted sum in their equn (32) corresponds to an average 
value of 460kms _1 . Clearly the more massive groups, in the 
tail of the distribution, dominate the deflection effect. If the 
same effect applies for the groups of Myers et al, the rms 
dispersion would rise to Srj = 0.45 in the coherent regime, 
as shown by the dotted line marked 'LM groups (Myers)' in 
Fig. 1. 

This is the extra dispersion in the CMB anisotropy 
caused by coherent scattering. Following Lieu & Mittaz we 
calculate the scale above which the scattering of light by the 
galaxy groups becomes incoherent rather than coherent. In 
the EdS model the angular diameter distance at the average 
group distance of z — 0.25 is cLa = 500h _1 Mpc. Therefore 
our assumed group diameter of 4h _1 Mpc subtends an angle 
of ss 0.45deg at this redshift. At smaller 9, Srj is constant, 
and at larger angles <5n oc 1/8 (see Fig. 1). Without the effect 
of the tail, the results are similar to the fl m = 1 results of 
Seljak (1996) although the <5n extends to larger scales before 
starting to decrease. With the tail, the model gives consider- 
ably higher values for <5n (see Fig. 1, dotted line) than found 
by Seljak for the standard ACDM model (long-dashed line). 

4.2 Galaxy masses from QSO lensing applied to 
results of Seljak (1996). 

We next apply the galaxy lensing results of Myers et al 
(2005) to the results of Seljak (1996). For the standard A cos- 
mology, Myers et al found values of the bias b m 0.05—0.15 at 
sub-Mpc scales, depending on the methodology. In the EdS 
cosmology, they found values in the range b w 0.15 — 0.3. 
Little evidence of scale dependence was seen in the data, so 
these results may be assumed to also apply at larger scales. 
Scaling a/8 from the unbiased results of CMBFAST in the 
A and EdS cosmologies according to the above derived anti- 
bias, the results for a/8 clearly rise by 1/6 i.e. by w lOx 
in the A case and ~ 5x in the EdS case. Clearly these re- 
sults are only approximate. In a future paper we shall relate 
the QSO galaxy lensing result of Myers et al (2005) to the 
CMB result more directly via equations (A 14) of Myers et 
al (2005) and (7) of Seljak (1996). 
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Figure 2. The C; for the O m = 1, h=0.35 baryonic model, as- 
suming reionisation with optical depth, r = 0.4 (dashed line). The 
lensed version of the C; with r = 0.4 is also shown (solid line) 
assuming the constant a (=0.005 rad.) model shown in Fig. 1. 
The simple SZ model described in Section 5 is also shown (long- 
dashed line). The data points are the 3-year WMAP results of 
Hinshaw et al (2006). 

Although the amplitudes have risen the form of the lens- 
ing dispersions remains the same. In particular, they flatten 
at small 8 much faster than the constant a model postulated 
in Section 3 above as having the appropriate form and am- 
plitude required to 'move' the first acoustic peak. The a/8 
result for £(r) = (r/7.4/i -1 Mpc) -3 is shown in Fig. 1 where 
the amplitude has been chosen to match the SCDM model 
at large scales. It can be seen that to reach the a = constant 
model given in Fig. 1, the r -3 mass correlation function may 
even have to steepen further to ensure a/8 maintains its — 1 
slope down to 7' scales, corresponding to ~ Hi -1 Mpc at an 
average galaxy depth of z ~ 0.15. The standard CDM model 
starts to flatten earlier at about 1 degree and the ACDM 
model at yet larger scales. This means that perhaps only a 
baryonic model with dissipation may produce the relatively 
steep mass correlation function in the r < 10h _1 Mpc range 
of interest for modifying the first acoustic CMB peak via 
lensing. 



5 APPLICATION TO THE MODEL OF 
SHANKS (1985) 

An example of a model which initially has the first acoustic 
peak at I — 330 rather than I = 220 is the Einstein-de Sitter, 
baryon dominated, low Ho model of Shanks (1985); Shanks 
et al. (1991); Shanks (2005). The C; from this model also 
produces too high a peak height, so to obtain a fit to the data 
we need to assume that the optical depth to reionisation is 
t — 0.4 which is at the high end of what has been suggested 
by the large-scale WMAP polarisation results (Kogut et al 
2003, Page et al 2006). In fact as shown by CMBFAST, 
as the t goes up, the polarisation power spectra quickly 
saturate and it becomes more difficult to reject higher values 
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of t. The resulting model Ci is compared to the WMAP data 
in Fig. 2. 

We then use the a = constant model from Fig. 1. There 
we have seen that this model has a high amplitude but one 
which may be motivated by the QSO magnification results of 
Myers et al (2003, 2005). The exponentially sharp increase in 
smoothing as I increases that then results from equation A7 
of Seljak is the sort of behavour which might be required to 
smooth the baryon model into the WMAP data. The result 
of applying this lensing smoothing function is shown as the 
solid line in Fig. 2 and it can be seen to give an improved 
fit to the first acoustic peak. (Note that we are assuming 
that l 2 C g i^{d) « 1 in equations A5 and A6 of Seljak in 
applying his equations A6 and A7 in this a ^constant case 
where a(6)/6 w 1 at 6 < 20'.) The question of whether the 
rise in a towards the smallest scales is self-consistent with 
the baryon model also arises. Certainly the linear theory 
baryon model does not produce such a rise, indeed showing 
less small scale power than the standard model. However, 
the non-linear power spectrum of the baryon model after 
top-down collapse is not easy to calculate and this leaves 
open the possibility that the baryon mass power spectrum 
may still be as steep as required to produce the a = constant 
model. As we have seen, a mass correlation function as steep 
as r~ 3 approximates the behaviour needed, compared to the 
observed galaxy correlation function's r -1 ' 8 behaviour. How- 
ever, galaxy correlation functions as steep as r~ 2 ' 2 are seen 
in samples of early-type galaxies on r < lh _1 Mpc scales (Ze- 
havi et al 2005, Ross et al 2006). The amplitude would have 
to be of order ro « 25h _1 Mpc implying strong anti-bias of 
the type suggested by the QSO magnification results. 

Note that we have only postulated the a = constant 
model in Fig. 1 and not derived it from equation (7) of Sel- 
jak (1996). Therefore, although the model has a/6 > 1 for 
6 < 20 arcmin, we have not violated the condition that 
equation (7) only applies in the case a/9 < 1 because of its 
assumption of the Born approximation. The main assump- 
tions we therefore have made in deriving the lensed Ci in 
Fig. 2 is that l 2 C g i^(6) << 1 in the application of Seljak's 
equation A7 as noted above. Nevertheless, the lensed result 
in Fig. 2 should still be regarded as only an approximation. 

It has also not been possible to fit the second peak. 
To address this issue, we could consider using the SZ ef- 
fect caused by the known presence of hot gas in the rich 
Abell clusters to re-create the second peak after its era- 
sure by the lensing. Myers at al (2004) have statistically 
detected SZ decrements around Abell clusters in WMAP 
data out to radii of 30 arcmin. These decrements simply 
add power to the power spectrum on cluster-sized scales. 
Assuming a random distribution of circular 'clusters', with 
a sky density of 3deg~ 2 , radius 12 arcmin and with 'decre- 
ment' AT/T = — 3 x 10~ 5 , we can produce a quite significant 
'second peak' (see Fig. 1). Although the match to the sec- 
ond peak is only approximate, this simple model with just a 
'top-hat' SZ profile at least serves to indicate the principle of 
this approach. The observed spectral index of the WMAP 
CMB fluctuations which effectively constrain the SZ con- 
tribution to the first acoustic peak (Huffenberger & Seljak 
2004) are much less effective constraints against a significant 
SZ contribution to the second peak. 

Further work is needed to check if the CMB EE and TE 
polarisation power spectra from WMAP can be reproduced 



by the above model (Kogut et al 2003, Page et al 2006). For 
example, it is known that lensing will produce B modes as 
well as E modes in the polarisation signal (Zaldarriaga & 
Seljak 1998) and it will be interesting to see if this provides 
a constraint on the current approach. 

For the moment, we conclude that using high-redshift 
reionisation to reduce the first two acoustic peaks, then lens- 
ing to move the first peak to larger scales and finally SZ 
to recreate the second peak, it appears possible to achieve 
an approximate fit to the WMAP CMB fluctuation power 
spectrum. Although this approach may be criticised as fine- 
tuning CMB foreground effects to achieve a fit to the data, 
it may also be viewed as being much less fine-tuned than 
the approach taken via exotic particles and dark energy, as 
employed to construct the standard ACDM model. 



6 CONCLUSIONS 

We first compared models for the effect of foreground galaxy 
clusters and groups on the CMB power-spectrum. We have 
shown that the lensing effects considered by Lieu & Mit- 
taz may have bigger effects than they suggested, if the full 
3-D velocity dispersion is used rather than the 1-D disper- 
sions assumed by these authors. This change produces « 3x 
bigger results for 6r)(= a(6)/6) than they suggested, which 
will increase the effect on the CMB peaks. 

We find that the empirical model of Lieu & Mittaz then 
gives significantly higher results for the magnification dis- 
persion at small scales (6 < 30') than those of Seljak (1996) 
assuming standard cosmological parameters. We have also 
shown that under the same foreground assumptions, now 
modelled in mass correlation function terms, the lensing for- 
malism of Lieu & Mittaz produces approximately the same 
results as the that of Seljak (1996). But both formalisms in 
the case of the standard ACDM foreground or in the em- 
pirical cases considered by Lieu & Mittaz tend to produce 
a/0 relations which are too flat towards small 6 and have 
too small amplitudes to modify significantly the first peak 
of the CMB power spectrum. 

In conjunction with a model that assumed lensing dis- 
persion a(9) = constant we then used the QSO magnifi- 
cation results of Myers et al (2003, 2005) to show that, if 
galaxy groups have bigger than expected masses, there could 
be significant effects of gravitational lensing on even the first 
peak of the CMB. Although the results of Myers et al remain 
controversial (Scranton et al. 2005, Mountrichas & Shanks 
2007), it is clear that at least in principle, lensing can af- 
fect the interpretation of the CMB power spectrum. But to 
have a significant effect on the first peak, the mass clustering 
power spectrum would then need to have been significantly 
underestimated in previous studies with anti-bias, b « 0.2, 
being required. 

There are many pieces of evidence which favour the 
standard model bias value of b « 1. For example, redshift 
space distortion analyses imply f3 — fiJ, 6 /6 w 0.5 ±0.1 
(e.g. Ratcliffe et al 1996, Hawkins et al 2003.) But such 
analyses tend to assume that the galaxies trace the mass 
with linear bias and in the case where groups contain as 
much mass as clusters this assumption does not apply. Even 
here we note that the infall velocity patterns around the 
rich environments of passive galaxies in 2dFGRS are very 
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similar to those around the poorer environments of active 
galaxies (Madgwick et al. 2003), suggestive of the evidence 
from the QSOs that groups dominated by spirals may have 
as much associated mass as the richer environments of the 
early- types. 

But still, the most direct previous constraints on the 
mass power spectrum come from estimates of cosmic shear. 
The QSO lensing results are again in clear contradiction 
with these since they are mainly consistent with b — 1. For 
example, Seljak (1996, Fig. 1) suggests that the cosmic shear 
upper limits on a(0)/9 are close to the SCDM predictions 
and well below the constant a model in Fig. 1. Yet even in 
the case of cosmic shear there are possible systematic ef- 
fects which may cause underestimates of the mass density 
and overestimates of b. For example, Hirata & Seljak (2004) 
have suggested that intrinsic galaxy alignments may prove 
difficult to disentangle from the effects of gravity on the 
shapes of galaxies. They suggest that even with the avail- 
ability of perfect redshift information to remove physically 
associated pairs of galaxies, the effect of lensing may still be 
corrupted by physical alignment if the foreground tidal field 
dynamically elongates the foreground galaxy while elongat- 
ing the background galaxy via its lensing effect. This could 
cause the cosmic shear to be significantly underestimated in 
present samples and explain the discrepancy with the QSO 
lensing result (see also Mandelbaum et al., 2006). 

Although the agreement in form between the galaxy and 
QSO power spectra and WMAP are impressive, there is the 
possibility that the galaxies do not trace the mass on the 
largest scales and may be affected by scale-dependent bias. 
The high masses for groups suggested by the QSO lensing 
results are certainly consistent with this idea. We also must 
recall that there are other measures of clustering which give 
much higher amplitudes of clustering than those given by the 
galaxies. For example, the rich Abell galaxy clusters have a 
correlation length of ro ~ 25h _1 Mpc. Another measure of 
clustering of mass at large scales is our large w 700kms _1 
motion with respect to the CMB (eg Kogut et al 1993). 
This has always been at the upper end of what is allowed 
if P « 0.5 and these results are made even more remark- 
able, given the observational evidence that the local region 
of ~40h _1 Mpc radius or more is also moving with a compa- 
rable 'bulk motion' (Lynden-Bell et al. 1988). 

Since the non-linear evolution of CDM models is well- 
worked through via collisionless N-body simulations, it is 
unlikely that such models can produce anything but the 
standard forms for a(9)/9. Dissipation in galaxy and cluster 
formation in the pure baryonic models could still leave the 
mass power spectrum steeper than expected from its initial, 
linear form. 

We have taken the view here that the crucial observa- 
tion is the high fraction of baryons in Abell clusters such as 
Coma. This 'baryon catastrophe' already forces the CDM 
model to include A to allow a model with a lower Q. m to be 
compatible with zero spatial curvature and inflation (White 
et al 1993). If the X-ray gas is to explain the missing mass 
in Coma, routes must be found to reconcile WMAP with 
the adiabatic Qbaryon = 1 prediction and we have shown 
that cosmic CMB foregrounds may form such a route. If the 
model can be shown to satisfy other constraints such as the 
WMAP polarisation power spectra, the result could be a 
model which is less finely-tuned than the standard model 



and is based on known physics, without recourse to unde- 
tected neutralinos or postulated dark energy. 
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